Fermi surfaces in general co-dimension and a new controlled non-trivial fixed point 
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Traditionally Fermi surfaces for problems in d spatial dimensions have dimensionality d — 1, i.e., 
codimension dc = 1 along which energy varies. Situations with dc > 1 arise when the gapless 
fermionic excitations live at isolated nodal points or lines. For d^ > 1 weak short range interactions 
are irrelevant at the non-interacting fixed point. Increasing interaction strength can lead to phase 
transitions out of this Fermi liquid. We illustrate this by studying the transition to superconductivity 
in a controlled e expansion near dc — 1. The resulting non-trivial fixed point is shown to describe 
a scale invariant theory that lives in effective space-time dimension D = dc + 1. Remarkably, the 
results can be reproduced by the more familiar Hertz-Millis action for the bosonic superconducting 
order parameter even though it lives in different space-time dimensions. 

PACS numbers: 



The ground state of a fluid of non-interacting fermions 
in spatial dimension d > 2 is typically a Landau fernii 
liquid. This state has a sharp Fernii surface and well de- 
flncd quasiparticle excitations. The Fernii liquid state 
can be understood within a fermionic renormalization 
group frameworkii^ in terms of a fixed point obtained as 
one focuses on modes within a bandwidth A of the Fernii 
surface and systematically reduces A. Phase transitions 
out of the Fermi liquid phase (such as the Stoner fer- 
romagnetic transition) are not described by traditional 
flows within this approach. An alternate approach of 
Hertz'', extended by Millis^, performs an RG analysis 
on the bosonic action for the order parameter obtained 
by integrating out the fermions. This bosonic action is 
nonanalytic in frequency and momentum and yields re- 
sults in agreement with those of Moriya^ who used a self- 
consistent renormalization approach. Since the approach 
of integrating out the gapless fermionic modes is contrary 
to the dogma of the RG that singularities should appear, 
not in the action but in the Green's functions, it has been 
a source of concern to some. 

In situations where the non-interacting band structure 
consists of distinct Fermi points (such as in graphene) 
the situation is much better. The low energy theory 
of such fermion systems is a massless 'relativistic' Dirac 
theory. In dimension d > 2 weak short ranged interac- 
tions are irrelevant in such a theory. At finite interaction 
strength various phase transitions that gap out the mass- 
less fermions can occur. Without the complications of an 
extended Fermi surface these can be analysed within the 
conventional field theoretic framework for critical phe- 
nomena. 

In this paper we consider a general class of problems 
that fall in between these two cases. The low energy the- 
ory of a system of non-interacting fermions in d spatial 
dimension may be characterized by the co-dimension of 
the surface in momentum space where the energy gap 
vanishes. The ordinary Fermi surface has co-dimension 1 
while the case of Fermi points has co-dimension d. It is 
interesting to study the general case of co-dimension dc 
with a eye toward gaining insight into the all important 
case dc = 1. In dimension d = 3 the case d^ = 2 is the 



one new possibility. This corresponds to fermions with 
line nodes in three dimensions which arise in a number 
of different situations, most strikingly in several uncon- 
ventional three dimensional superconductors. 

We first use the fermionic RG methods of Ref. |l||j 
to argue that whenever dc > I all short range interac- 
tions are irrelevant. With increasing interaction strength 
phase transitions which gap out the fermions are again 
possible. For dc < d these occur in the presence of an 
extended gapless 'Fermi surface' and thus share the com- 
plications of the usual case dc = I- As a concrete exam- 
ple we study a pairing phase transition associated with 
superconductivity. We show that this can be accessed 
within the fermionic RG through a controlled e expan- 
sion in the co-dimension dc near dc — I. Remarkably a 
non-trivial scaling structure with effective space-time di- 
mensionality dc+1 is obtained. The transition can also be 
analysed within the Moriya-Hertz-Millis approach which 
works with a bosonic order parameter that lives in the 
full d space dimensions. Despite the difference in dimen- 
sionalities both approaches give identical results! Com- 
parison of the two approaches thus provides much insight 
into the nature of quantum criticality in fermionic sys- 
tems with an extended gapless Fermi surface. 

Consider then a 'generalized Fermi surface' of dimen- 
sion d — dc- Coordinates on this Fermi surface play the 
role of internal symmetry degrees of freedom and do not 
play a big role in the low energy limit for weak interac- 
tion strengths. We refer to them as isospin directions. 
We will assume a dispersion where the energy varies lin- 
early with momentum deviation away from this surface 
(as happens for usual Fermi surfaces and for line nodes 
in 3d). Let us call these directions along which energy 
varies linearly with momentum the Dirac directions. 

We begin in d = 2 when we have a Fernii line (one 
isospin direction, parameterized by the angle 6) and one 
radial Dirac direction measured by A:. In this world the 
momentum transfer q is two-dimensional. We are inter- 
ested in q because any bosonic order parameter bilinear 
in fermions will have q and uj at its arguments. Let us 
now increase the overall dimensions by e so that we have 
one isospin direction, 1 + e Dirac dimensions and 2 + e di- 



rections for q. Along with w we have a total of D = 2 + e 
dimensions for the (w,k) vector (on which the fermion 
fields depend ) and D — 3 + e dimensions for the (w,q) 
vector on which the order parameter depends. The case 
e = 1 corresponds to our target, the Fermi line in d = 3 
with two Dirac dimensions. 

The action for this theory (with Fermi velocity vp = 1) 
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is invariant under the RG transformation: 


(1) 


A -> s-'^A 


(2) 


Oj' — SbJ 


(3) 


k' = sk 


(4) 


i^{0jM) = s(3+e)/2^/(^,k,g) 


(5) 



The angle 6 that parameterizes the Fermi surface does 
not renormalize. Note that in an actual problem the 
Dirac line (Fermi surface) may not be a circle. The uni- 
versal answers calculated here as well as the kinematical 
implications in the limit A/ Kp —> are insensitive to 
this. 

Let us now add to this the BCS interaction. Recall 
from Refsiii^ that in the limit A/Kf -^ 0, the require- 
ment that all four lines lie in the bandwidth and obey mo- 
mentum conservation means the only possible four- Fermi 
interactions are those with either nearly forward scatter- 
ing described by a function m(0i, 62) (which will become 
u{6i —62) for a circular Fermi surface) and those with 
nearly zero incoming momentum described by vifii^O^) 
(which will become v(6i — O3) for a circular fermi sur- 
face). We say "nearly" in both cases because deviations 
of order A are kinematically allowed. These amplitude u 
will always be possible since if the two incoming lines lie 
in the bandwidth, (nearly ) forward scattering will en- 
sure the two final lines satisfy momentum conservation 
and lie in the bandwidth. The same logic works for v 
as long as we pick one incoming and one outgoing line in 
the bandwidth, for their opposites in the Cooper pair will 
satisfy momentum conservation and lie in the bandwidth 
as long as iJ(— k) = -E(k). 

By power counting the interactions u and v that sur- 
vive the restriction to the narrow bandwidth near the 
Fermi surface scale as 



(6) 



and similarly for u. Thus these are irrelevant for all co- 
dimension dc > 1, and the low energy physics is that of 
the free fermion theory. 

Let us now study a phase transition that occurs at fi- 
nite interaction strength that gaps out all the fcrmions. 
We consider a superconducting transition driven by in- 
creasing V. We consider a simple model where v is con- 
stant except for the antisymmetry demanded by the Pauli 
principle: 
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where labels like (2) are shorthand for angles, momenta 
and frequency labeled 2, and there is no integral over (4) 
which is required to equal (1 + 2 — 3) by conservation laws. 
If V is large the fermions will pair and a superconducting 
state will result. We now study the transition to this 
state to leading order in an e expansion for small e. 

If we eliminate the modes between A and A/s we get 
an all too familiar loop correction 
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where the loop may be done in D = 2 (which is e = 0) 
since we expect the v in front of it to be of order e at the 
fixed point. Thus 
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We can now read off the correlation length exponent: 

„=l/(3'{v*)^l/e. (11) 

We need to find the exponent that corresponds to the 
"magnetic field" that couples to the order parameter 
A(a;, q) — ((V'V'))c^.q- So let us consider a source term 
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where —9 is the angle on the Fermi surface where the 
momentum is opposite to that at the point 9. 

In this expression, one field Tp lies inside the bandwidth 
and the other has the opposite frequency and momentum 
plus a tiny w, q. It is readily verified that at tree level 



h' = hs 
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independent of e. At one loop, the source h couples to tptp 
via a particle-particle bubble (Figure 1) whose momenta 
lie in the shell being eliminated. Doing the calculation at 
e = due to the v* that multiplies the bubble, we find 



h' = hs{l + v*\ns) 
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All we need now is the scaling law for E^mg the singu- 
lar part of the ground state energy per degree of freedom. 
Recall the argument. If we eliminate some high energy 
variables in a path integral, it will make an additive non- 
singular contribution the free energy, to be carried on 
to the next stage of mode elimination. Since the singu- 
lar part is the same, the singular part of the energy per 
degree of freedom will scale as 
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FIG. 1: Renormalization of the field that couples to the BCS 
order parameter at one loop. Open (closed) circle is the bare 
(renormalized) vertex . 



where Sv = v — v* and D is the number of dimensions 
along which we have thinned out degrees of freedom and 
rescaled frequency or momentum. The answer in our case 
is clearly D = 2 + e. Henceforth the subscript in Esmg 
will be dropped. Thus we write 
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By taking appropriate derivatives with respect to h and 
5v and then setting s = {Sv)~^^'^ = {Sv)~'^ , h = we find 



a 
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r] = 2-e 
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V = l/e 
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Some comments are in order. By a we mean the expo- 
nent for the second derivative of E with respect to the 
control parameter 5v used to tune the transition. Note 
that a < 0. The exponents /3 and 7 come from the first 
and second derivatives with respect to h. One can find 
•q, the anomalous exponent for the order parameter cor- 
relation by directly computing it to one loop at w* = e 
or by using the relation 
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We emphasize that the fermionic RG gives a canon- 
ical description of the critical point. In particular con- 
ventional hyperscaling equations are satisfied with space- 
time dimensionality dc + 1. 

Let us ask what we would find had we chosen the 
Hertz-Millis route of integrating out the fermions to get 
an action for the pairing field A. We would find 

5a = I A*(^,q)((<z2 +0,2)^/2 +r)A(c.,q)^d2+-q 

+uA*AA*A terms-l- ... (23) 

Note that q lives in 2 -I- e dimensions and we have a 
D — 'i + e theory. It has to reproduce a fixed point that 
we saw is of effective dimensionality D = 2 + e. Let us see 
how far we can go with the quadratic part of the action 
which is just the Cooper pair susceptibility. 

First, the tree level scaling for the critical r = theory 
is 

A = s(3+2^)/2A' which means h' ^ s^^+^^^'^h (24) 



from which it follows that 



r = rs 



which means v 
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in agrement with the old answer, even though A scales 
differently now. 

From Eq. [M] it follows that the quartic coupling u 
in Eqn. [23] is highly irrelevant ij it is nonsingular in 
its to and q dependence. But it is singular due to the 
fermion loop that generates it, but still irrelevant. Let 
us now calculate [3 paying attention to this dangerously 
irrelevant term. 

If we evaluate the fermionic loop contribution to u, 
with all external q and w equal to zero we will find (from 
a box diagram of 2 pairs of fermions of equal and opposite 
('^,k)). 



dwdi+'^k 



(tj2 + fc2)2 

The infrared divergence of this diagram tells us that 
M(cj,q) - 
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(27) 



(Vc^2+^2)2-e- 

From this and Eqn. [23] we may deduce that 

uo(s) = s-iuo- (28) 

As in the case of (f)^ theory above 4 dimensions the u 
will be dangerously irrelevant and it will be necessary to 
keep it to obtain correct results for some exponents. 

Turn now to the energy density which scales as follows 
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The exponents v, 7 and 77 can be deduced by working in 
the non-superconducting state {i.e r > 0). Then we may 
safely set u — 0, and the results are identical to that 
obtained above by the fermionic RG. 

The exponents /3 and a require retaining the irrelevant 
interactions. From the first derivative of E with respect 
to h we obtain the order parameter m —< A >: 
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and the mean-field result far from the transition when 
rs"^ is of order unity 
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Setting rs^ = 1 we find 



p^l/e 



(31) 



(32) 



again in agreement with the fermionic RG. Finally a (ex- 
ponent for the second derivative with respect to the con- 
trol parameter) may be obtained within mean field theory 
by examining the ground state energy in the presence of 
a uniform condensate A which assumes the singular form 
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in terms of r and u which are viewed as fixed (and not 
running) parameters. At the minimum A ~ (r/u)^'^ 
(so that (3 = 1/e) while the second r -derivative of the 
minimum energy gives a = 1 — - , in agreement with the 
fermionic RG. Note that the Gaussian action for A in 
the disordered phase describes the fluctuation correction 
which is subdominant to the mean field answer. 

The agreement of the fermionic and bosonic versions is 
remarkable and seems unlikely given that they describe 
theories in different number of effective dimensions (2 + e 
for fermions and 3 + e for bosons). Despite this and 
the fact that A scales differently in the two versions, all 
exponents are the same order e. It is always possible that 
something differs at higher orders in the s expansion. 

There are reasons for believing that the one-loop order 
e results are exact. As stressed in Ref.^ a narrow band- 
width problem has a small parameter A/Kp that plays 
the role of 1/A^. It is well known that in the large N limit 
the one- loop beta function is exact. Another way to say 
this is that one may simply do the RPA. As explained in 
Refi^ the RPA sum over repeated bubbles becomes exact 
because there is no phase space for any other diagram in 
the hmit A/Kp -^ 0. 

Despite this agreement we emphasize that the 
fermionic approach is much more natural as it properly 
gives a scale invariant fixed point with a correct effective 
dimensionality of dc + 1- In particular all hyperscaling 
equalities are satisfied with this effective dimension. The 
Hertz-Millis theory does not, at first sight, look like a 
theory that will satisfy any such scaling. However after 
including for the complications of the dangerously irrel- 
evant interactions this hidden effective dimensionality is 
correctly captured. Indeed the Hertz-Millis description 
needed to be more complicated in order for such a d -t- 1- 
dimensional theory to scale as though it lived in a differ- 
ent lower dimension dc + 1. 

Though we will not discuss it here clearly we can now 
go ahead and compute universal physics associated with 
the transition at finite temperatures to study the quan- 



tum critical region^. One aspect of this theory is worth 
considering. What is the fate of the fermions themselves 
right at the critical point? This may be straightforwardly 
studied by considering the self-energy of the fermions at 
the fixed point. We find that the fermionic quasiparti- 
cles are not destroyed even at this non-trivial interacting 
fixed point. In other words the quasiparticlc residue Z 
stays non-zero all the way upto and including the crit- 
ical point. The disappearance of the generalized Fermi 
surface on entering the superconducting state thus hap- 
pens without ever destroying the fermionic quasiparticlc. 
This is a different route to killing a Fermi surface from 
the Mott-like transitions considered in Ref. ,81 where Z 
vanishes at the critical point^. Note that the Fermi sur- 
face stays sharp at the present critical point as expected 
on general grounds^. 

To conclude we studied fermion systems with gener- 
alized 'Fermi surfaces' with co-dimension dc in d space 
dimension in the presence of interactions. Not surpris- 
ingly weak short range interactions were irrelevant for 
all dc > I. We studied a non-trivial superconducting 
transition out of such a phase through a controlled co- 
dimension expansion using the fermionic RG technique. 
The resulting theory scales with effective dimensionality 
dc + I. We also demonstrated that this result can be 
recovered within the Hertz-Millis approach despite the 
different aparent dimensionality. More generally study- 
ing phase transitions of Fermi surfaces with general co- 
dimension dc in d space dimensions is a potentially useful 
theoretical device. In the future it might be interesting to 
understand the Stoner transition or problems with gauge 
fields by generalizing to dc > 1. Perhaps such generaliza- 
tions will lead to controlled expansion methods suitable 
for such problems. 

We are grateful to Ganpathy Murthy for very use- 
ful conversations and the NSF for grants DMR-0705255 
(TS), and DMR- 0103639 (RS). RS is grateful to the Cen- 
ter for Condensed Matter Theory at MIT for hosting his 
sabbatical in Spring 2008, which made this work possible. 



^ The RG for fermions was briefly addressed in G. Benfatto 
and G. Gallavotti, Phys. Rev. B 42, 9967, (1990), J. Feld- 
man and E. Trubowitz, Helv. Phys. Acta 63 (1990) 157, 
R. Shankar, Physica A177, 530 (1991); J. Polchinski, , J. 
Polchinski, Effective Field Theory and the Fermi Surface, 



http://arxiv.o rg/abs/hep-th/9210046v2 
The RG for fermions was addressed in great detail in 
R.Shankar, Rev. Mod. Phys. 66, 129 (1994). 
^ T. Moriya, Spin Fluctuations in Itinerant Electron Mag- 
netism, Springer- Verlag, Berlin (1985); 



'' J. A .Hertz, Phys. Rev. B 14, 1165, (1976). 

^ A. J. Millis, Phys. Rev. B 48, 7183, (1993) 

® J. Polchinski, Low energy dynamics of the spinon-gauge sys- 
tem, cond-mat 9303037 

^ S. Sachdev, Quantum Ptiase Transitions, Cambridge Uni- 
versity Press, (1999). 

* T. Senthil, arX iv:0803.4009l 

^ W. F. Brinkman and T. M. Rice, Phys. Rev. B 2, 4302 
(1970). 



